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I  We  prove  from  quantum  electrodynamics  that  two-photon 

emission  in  transitions  between  atomic  inner  shells  can  be 
treated  in  lowest-order  approximation  in  the  single-particle 
model.  taking  many-r^body  effects  into  account  >f^show  that,  in 
the  reduced  single-particle  model,  the  Pauli  exclusion  principle 
does  not  prohibit  summing  over  all  possible  intermediate  states, 
including  core  states.  Some  transition-rate  formulas  are 


presented  that  are  useful  for  numerical  calculations 


AFOSR-TH.  R  7  -  1  42  4 


I.  INTRODUCTION 

Recent  developments  in  experimental  techniques  have  made  it 
possible  to  measure  two-photon  transitions  between  atomic  inner- 
shell  vacancy  states.^  The  results  must  be  interpreted  relativis- 
tically.  Existing  theoretical  work  has  been  based  on  the  hydroge- 
nic^”"^  or  nonrelativistic  self-consistent-field  (SCF)^  ^  models, 
except  for  the  recent  relativistic  SCF  calculation  of  Mu  and 

7  , 

Crasemann.  The  latter  work  suggests  exploration  of  some 
interesting  questions:  Is  the  Pauli  exclusion  principle  violated 
if  one  includes  occupied  orbitals  among  the  intermediate  states, 
as  done  in  Ref.  7?  What  are  the  theoretical  similarities  and 
differences  between  two-photon  processes  in  hydrogenic  vs.  inner- 
shell  atomic  systems?  In  this  paper  we  address  these  questions. 

In  Sec.  II,  we  analyze  the  theory  of  two-photon  transitions  on 
the  basis  of  quantum  electrodynamics  (QEO) .  The  possible 
intermediate  states  are  identified  automatically;  they  can  be 
classified  into  two  types  according  to  time  ordering.  One  type  is 


advanced:  the  electron  final  state  occurs  when  the  intermediate 


far 


states  do,  which  is  in  advance  of  the  disappearance  of  the  ini-  I  ^ 

□ 

tial  state.  The  second  type  is  retarded:  the  time  order  for  the  d  □  i 

Ion _ 

electron  is  first  initial  state,  then  intermediate  states,  then - — 


final  state.  In  this  derivation,  the  Pauli  exclusion  principle  ison/ 


built  in  by  commutation  relations  of  operators.  The  final  expres- 


itr  CodoB 


and/or 
.ipuclai 


sion  for  the  S-matrix  element  shows  that  the  problem  can  be 
formally  treated  as  a  time-independent  one,  with  the  electron 
moving  in  the  field  of  the  nucleus  and  the  other  electrons.  In 
the  reduced  model,  since  there  is  no  time-order  difference,  the 
original  two  types  of  intermediate  states  patch  up  a  new  set  of 
intermediate  states  which  is  a  complete  set.  The  Pauli  principle 
is  thus  seen  to  allow  summation  over  all  possible  intermediate 
states,  including  occupied  inner-shell  states. 

In  Sec.  Ill,  a  time- independent  many-body  approach  is 
delineated.  In  Sec.  IV,  as  an  application  of  the  results  of 
earlier  sections,  we  provide  formulas  for  matrix  elements, 
transition  rates,  and  differential  equations  to  be  solved,  which 
are  formally  the  same  as  for  the  hydrogenic  case.  The  only 
difference  is  the  external  potential  acting  on  the  electron. 

It  is  not  that  an  inner-shell  vacancy  is  trivially  equivalent  to 
the  hydrogenic  case:  there  are  many-body  effects  in  inner  shells, 
particularly  the  restrictions  imposed  by  the  Pauli  exclusion 
principle,  while  the  hydrogenic  system  consists  of  just  a  single 
bound  electron.  Nor  should  the  reader  think  that  this  paper's 


result  is  just  a  simple  repetition  of  the  conclusions  of  Brown  et 

O 

al.  regarding  gamma-ray  scattering,  which  are  based  on  a 
single-electron,  i.e.  hydrogenic,  model.  If  we  consider  many-body 
effects  and  the  built-in  Pauli  exclusion  principle  in  inner-shell 
gamma-ray  elastic  scattering,  the  core  contributions  should  in 
fact  be  excluded  when  summing  over  intermediate  states  (unlike 
the  procedure  in  Ref .  8) .  By  contrast,  as  shown  here,  all 
intermediate  states  must  be  included  (as  done  in  Ref.  7)  when  one 
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calculates  two-photon  inner-shell  transition  probabilities, 
taking  the  Pauli  exclusion  principle  into  account. 


II.  QUANTUM-ELECTRODYNAMIC  APPROACH  AND  THE  QUESTION  OF 

INTERMEDIATE  STATES 

The  dominant  part  of  two-photon  emission  is  of  the  double 
electric-dipole  (2E1)  type;  this  is  the  only  kind  considered 
here. 

In  the  theory  of  QED,  the  second-order  scattering  operator  S 
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can  be  expressed  in  Furry's  representation  as 


The  matrix  element  is 

<o|<corelc.  (2) 


where  |core>,  in  the  lowest-order  approximation,  stands  for  an 
antisymmetrically  filled  core  wave  function,  and  |0>,  for  a  photon 
vacuum  state.  We  denote  photon  creation  and  annihilation 
operators  by  c^  and  c,  respectively,  and  electron  creation  and 
annihilation  operators  by  a^  and  a.  The  remainder  of  the 
notation  is  of  standard  usage.  Natural  units  (tl=c=l)  are 
employed  throughout.  OC  =  i/f}7.  0S6 . 

The  electron  field  operator  can  be  expressed  as 


(3) 


The  summation  index  a  runs  over  all  core  states,  and  n,  over  all 
the  excited  states. 


After  time-ordered  contraction  manipulation,  the  S  matrix 
exement  takes  the  form 

5^,  =  -  tx;  f Am  A',*  t  aI;,  a',*) 


where 


x')=z  -i  <corcl  T  ^  1  corc> 


is  the  propagator  for  an  electron  hole,  andl(^(x)  and^(x)  are 
bispinor  wave  functions.  The  indices  1  and  2  denote  the  initial 
and  final  states  of  electrons,  respectively. 

The  explicit  form  for  the  electron  hole  propagator  can  be 
shown  to  be 


^  ,,fL  £-£,  +  n 


The  first  term  in  the  bracket  on  the  right-hand  side  of  Eg. 
(6),  which  involves  core  intermediate  states,  is  an  advanced 
Green's  function  which  describes  the  t  >  t'  process;  the  second 
term  which  involves  excited  intermediate  states,  is  a  retarded 
Green's  function  which  describes  the  t'  >  t  process.  This 


distinction  is  important  in  understanding  the  physical  meaning  of 
the  two  types  of  intermediate  states. 

The  two-photon  processes  can  be  described  by  Feynman 
diagrams.  We  specify  the  initial  and  final  hole  states  in  the 
following  examples.  Time  progresses  from  left  to  right.  A  solid 
line  with  arrow  opposite  to  the  time  direction  represents  an 
electron  hole,  a  line  with  arrow  in  the  time  direction  represents 
a  real  electron.  There  are  two  eguivalent  ways  of  expressing  the 
advanced  process,  as  shown  in  Fig.  1.  Diagram  1(a)  describes  the 
process  in  terms  of  a  hole,  and  diagram  1(b),  in  terms  of  an 
electron.  For  the  retarded  process,  there  are  two  similar 
diagrams,  shown  in  Fig.  2. 

From  this  analysis,  we  obtain  only  two  types  of  intermediate 
states.  Two  other  types  of  intermediate  states  could  in 
principle  be  constructed,  but  do  not  occur  in  our  derivation.  One 
of  these  pertains  to  retarded  real  core  states;  the  disappearance 
of  this  type  embodies  the  Pauli  exclusion  principle.  The  other 
type  of  intermediates  that  does  not  appear  here  involves  advanced 
excited  hole  states;  this  type  is  excluded  by  the  fact  that  an 
excited  state  cannot  be  annihilated  from  the  core. 

Electromagnetic  waves  and  electron  waves  have  the  forms 


-iw.t 


(7) 


,  v;(x)= ’h(>:) 


respectively.  After  we  integrate  over  time  t  ,  then  t',  and  then 


energy  ,  we  find  the  following  expression  for  the  S  matrix 


element: 


-i  2ri«ijd>x  fl')  >” 


a  w,-£.t£,»i1  n  f.+u,-E.*>l  (8) 

f  w.-£.t£4*i-l  •  J 


The  function  factor  expresses  the  energy  conservation  law 

(9) 

By  using  this  energy  conservation  condition,  we  can  rewrite  the 
energy  denominators.  It  then  becomes  possible  to  combine  the  two 
types  of  Green's  function  into  one  with  a  complete  set  of  wave 

functions  as  its  intermediate  states; 

5.,=  -iZRotJJdJr 

^  (10, 

[i  f  £,  +  u).-E.  J 

<rfe, -w*) 

Here,  the  index  i  runs  over  all  possible  intermediate  states. 

The  expression  for  the  S  matrix  element  is  thus  reduced  to  a 
time-independent  single-particle  version,  which  is  formally  the 
same  as  in  the  hydrogenic  model  if  we  call  the  electron  2  as  the 
initial  state,  the  electron  1  as  the  final  state. 


The  only  differences  between  the  two  models  are  in  two 
aspects:  (1)  the  hydrogenic  electron  moves  in  the  Coulomb  field 

of  the  nucleus,  while  the  electron  displaced  by  the  inner-shell 
hole  moves  in  the  combined  field  of  the  nucleus  and  the  remaining 
electrons.  (2)  The  summing  and  averaging  over  initial  and  final 
states,  respectively,  are  opposite  in  the  inner-shell  case  to 
that  of  outer  shells.  We  have  thus  proved  a  theorem  which  is 
useful  in  deriving  formulas  for  numerical  calculations: 

THEOREM  I. 

Two-El  photon  emission  from  atomic  inner  shells,  filling  one 
vacancy,  can  be  expressed  in  the  lowest-order  approximation  in  a 
single-particle  formulation  which  is  formally  the  same  as  in  the 
hydrogenic  case,  with  the  electron  hole  being  treated  as  a  real 
electron.  The  only  differences  are  in  the  potential  acting  on  the 
electron  and  the  treatment  of  summing  and  averaging  over  the 
initial  and  final  states. 

III.  TIME-INDEPENDENT  MANY-BODY  APPROACH 

In  a  time-independent  many-body  approach,  the  initial  and 
final  states  can  be  written  as 

|i  >=  a.,  lcorc>  ,  |2>  =  d,l  core;>>  (H) 

As  discussed  in  Sec.  II,  we  have  two  types  of  intermediate 
states,  advanced  and  retarded.  In  time-independent  many-body 
perturbation  theory,  these  can  be  written 


a*  I  core> 


(12) 


aX  a.,  a,  | 


The  dipole  operator  has  the  form 


S"  =  £  9m  a*  3ij 
■j 

where  m  is  the  magnetic  quantum  number.  The  time-independent 
many-body  expression  for  the  emission  rate  is 

‘  E,-£v-<o.  i 


After  contractions,  this  form  is  reduced  to  a  single-body  form 


.nij  a"’*  a  .m, 

ciw  1  I  jix*  T-  ^  r  ***  c-  **» 

—  =  fif  *  r-t,tet-u^  «  £.-£,-“<  §  £.-£«-ai. 


*•» 


,a 


The  many-body  forms  of  the  two  types  of  intermediate  states 
are  now  reduced  to  core  states  and  excitations  in  the  single-body 
form.  If  we  impose  the  energy  conservation  law  Eg.  (9) ,  the 
Green's  function  that  involves  core  states  and  the  one  that 
involves  excitations  are  combined  into  a  new  Green's  function 
that  involves  a  complete  set  of  intermediate  states.  This  proves 
Theorem  I  over  again.  We  therefore  have  the  final  form  for  the 


emission  rate: 


d\v  1  j  y  y  ^‘*  +  T  — — 

—  =  w,  to.  1 


If  we  specify  the  initial  state  for  the  atom  as  a  [Is]  hole  state, 
the  possible  final  states  will  be  [2s],  [3s],  [3d  ],  [3d  ], 

J/ ^  D/ ^ 

[4s],  .  .  .  ,  hole  states.  The  possible  intermediate  states  are 

of  two  types.  The  ones  derived  from  retardation  are  of  the  type 

of  single-hole  p  states  of  the  core,  and  the  ones  derived  from 

advancement  are  of  the  type  of  the  states  with  one  excited  p 

electron  and  two  holes  in  the  core.  This  agrees  with  the 
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classification  of  Aberg,  and  of  Freund  and  Bannett.  E.J. 
Robinson,  M.  Crance  and  M.  Aymar^^  had  similar  treatments  and 
conclusions  for  core  excitations  in  outer  shell  problems. 

The  time- independent  many-body  approach  leads  to  the  same 
result  as  the  QED  approach.  It  is  simpler,  but  involves  some 
imposed  conditions,  such  as  those  regarding  the  intermediate 
states  and  the  energy  conservation  law.  It  cannot  reveal  the 
time-ordered  process,  and  hence  leads  to  less  physical  insight. 

IV.  MATRIX  ELEMENTS  AND  EMISSION  RATES 
Theorem  I  reduces  the  inner-shell  vacancy  two-photon  decay 
problem  to  one  tractable  in  the  single-particle  model.  The 
formulas  for  the  hydrogenic  case  can  all  be  adapted  to  the  many- 
body  problem,  but  the  electron  wave  functions  must  be  calculated 
for  different  potentials.  A  formalism  for  relativistic  two- 
photon  emission  in  the  hydrogenic  case  has  been  given  by  Goldman 


dW 


—  27Ii  d-  y  E-J 


The  transition  amplitude  is  given  by 


*  '^0 

where  S^^»and  T^^iare  solutions  of  the  perturbed  Dirac  equation 

(m-S,  ^cj,  *V(r>)S^^.(r,  w.)  +  (r,uj,) 

.  ,  < 


a 


The  inhomogeneous  terms  K^^.and  L^^>and  the  terms  and  in  the 

integrand  of  Eq.  (22)  are 

+  (coulomb  gauge) 

.  I  <*  Wi  r  f*j%T 


(r,w.)  =  [  ^  *^>r  ''  «.r 


^  fi  (r)  •*■  r -  l) 


(length  gauge) 


(Coulomb  gauge) 
(length  gauge) 


(Coulomb  gauge) 
(length  gauge) 


(Coulomb  gauge) 
(length  gauge) 


Here,  G^(r)  and  (r)  are  the  large  and  small  components  of  the 
electron  wave  function,  satisfying  the  radial  equation 


(Tn-  -  ■T)F^(r)=0 

-  f ^  (O  -  VM)F^(r)  =  0 


These  formulas  are  formally  identical  to  those  for  the  hydrogenic 

case.  The  only  difference  is  in  the  potential  V(r) .  In  Eqs. 

(24)  and  (25) ,  V(r)  is  the  potential  due  to  the  nucleus  and  the 

other  electrons,  but  we  do  not  exclude  the  core  contribution  in 

the  inhomogeneous  equations  (22).  In  the  hydrogenic  case,  V(r) 

is  the  pure  Coulomb  potential  due  to  the  nucleus.  The  notation 

.  2 

used  here  is  the  same  as  that  in  Parpia  and  Johnson's  paper. 

If  we  set  K'  =  -1  in  Eg.  (23),  we  obtain  the  s-*s  decay  rate,  which 

2 

is  the  same  as  the  formula  of  Parpia  and  Johnson.  If  we  set 
K’=  2,  and  -3,  we  find  the  ds/j-^s  and  dy/a-«  decay  rates. 


ACKNOWLEDGMENTS 


This  work  was  encouraged  by  W.  R.  Johnson.  The  author 
thanks  B.  Crasemann  for  stimulating  discussions  and  helpful 
comments  on  the  manuscript.  The  author  acknowledges  fruitful 
discussions  with  Xiao-Ling  Liang,  Xing-Dong  Mu  and  T.  Aberg. 
This  work  was  partially  supported  by  the  Air  Force  Office  of 
Scientific  Research  under  Grant  87-0026  and  by  the  National 
Science  Foundation  under  grant  PHY-85  16788. 


K.  Ilakovac,  J.  Tudor ic-Ghemo,  B.  Busic  and  V.  Horvat,  Phys. 
Rev.  Lett.  56,  2469  (1986). 

W.  R.  Johnson,  Phys.  Rev.  Lett.  29,1123  (1972);  F.  A.  Parpia 
and  W.  R.  Johnson,  Phys.  Rev.  A  26,  1142  (1982) . 

S.  P.  Goldman  and  G.  W  .F.  Drake,  Phys.  Rev.  A  24,  183  (1981) 
V.  Florescu,  Phys.  Rev.  A  30,  2441  (1984) . 

T.  Aberg,  in  Atomic  Inner  Shell  Processes,  edited  by  B. 
Crasemann  (Academic,  New  York,  1975),  p.353. 

I.  Freund,  Phys.  Rev.  A  7,  1849  (1973);  Y.  B.  Bannett  and  I. 
Freund,  Phys.  Rev.  A  30,  299  (1984). 


7.  X.  Mu  and  B.  Crasemann,  Phys.  Rev.  Lett.  57,  3039  (1986). 

8.  G.  E.  Brown,  R.  E.  Peierls,  and  J.  B.  Woodward,  Proc.  R. 
Soc.,  Ser.  A  227,  51  (1954). 

9.  V.  B.  Berestetskii,  E.  M.  Lifshitz  and  L.  P.  Pitaevskii, 

Quantum  Electrodynamics ,  (2nd  Ed.,  p.  482.,  1982). 

10.  E.J.  Robinson,  Phys.  Lett.  A  117,  465  (1986). 

11.  M.  Crance  and  Mireille  Aymar  J.  Phys.  B20.  L155  (1987). 


Feynman  diagrams  for  advanced  process 

(a)  in  terms  of  a  hole,  (b)  in  terms  of  an  electron 

Feynman  diagrams  for  retarded  process 

(a)  in  terms  of  a  hole,  (b)  in  terms  of  an  electron 


